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A computational conjugate heat transfer methodology was developed, as a first step 
towards an efficient and accurate multiphysics, thermo-fluid computational methodology to 
predict environments for hypothetical solid-core, nuclear thermal engine thrust chamber 
and components. A solid conduction heat transfer procedure was implemented onto a 
pressure-based, multidimensional, finite-volume, turbulent, chemically reacting, thermally 
radiating, and unstructured grid computational fluid dynamics formulation. The conjugate 
heat transfer of a cylindrical material specimen heated by an impinging hot hydrogen jet 
inside an enclosed test fixture was simulated and analyzed. The solid conduction heat 
transfer procedure was anchored with a standard solid heat transfer code. Transient 
analyses were then performed with variable thermal conductivities representing three 
composites of a material utilized as flow element in a legacy engine test. It was found that 
material thermal conductivity strongly influences the transient heat conduction 
characteristics. In addition, it was observed that high thermal gradient occur inside the 
cylindrical specimen during an impulsive or a 10 s ramp start sequence, but not during 
steady-state operations. 


Nomenclature 

C h C 2 ,C 3 ,C^ turbulence modeling constants, 1.15, 1.9, 0.25, and 0.09. 
C p = heat capacity 

D = diffusivity 

H = total enthalpy 

K = thermal conductivity 

k = turbulent kinetic energy 

p = pressure 

Q = heat flux or heat source 

T = temperature 

t = time, s 

Uj = mean velocities in three directions 

x = Cartesian coordinates 

e = turbulent kinetic energy dissipation rate 

0 - energy dissipation contribution 

p = viscosity 
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Hi = turbulent eddy viscosity (=pC^k 2 /s) 

T1 = turbulent kinetic energy production 

p = density 

a = turbulence modeling constants 

r = shear stress 

co = chemical species production rate 

Subscripts and superscripts 
cl = centerline 

s = surface 

t = turbulent flow 

v = volume 

0 = represents cells at the fluid/solid interface 

1 = represents cells at one point into the solid 

= represents the next time-level solution 


I. Introduction 

N uclear thermal propulsion (NTP) may open up the solar system to far broader and faster exploration than is 
now possible with chemical propulsion. The feasibility of solid-core NTP systems was established by extensive 
testing in the Rover/NERVA programs and the technical merits of NTP have been identified and summarized. 1 The 
basic concept involves a reactor core consisting of several hundred flow elements that are used to heat a low 
molecular weight gas, such as hydrogen. Thrust is generated as the hot hydrogen expands through a high expansion 
ratio nozzle. The need to push flow element temperature to extremes in order to maximize performance intensifies 
the possibility of material failures. 2 For example, the mid-section corrosion was named after the materials at the 
mid-section of the flow element was found cracked after tests. 2 In order to develop candidate high temperature 
materials that would be compatible with the hot-hydrogen environment of a high performance solid core NTP 
engine, new materials for flow elements need to be researched and tested. And accompanying computational 
methodology needs to be developed for analysis of test results in order to understand the contributing factors that 
lead to cracked materials, and to develop steps to improve the materials and operating procedures, so as to avoid 
future component failures. 

In this work, a computational conjugate heat transfer methodology was developed to study the effect of thermal 
conductivity on thermal gradient inside a cylindrical specimen in a hot hydrogen environment, paralleling a hot 
hydrogen materials test effort that will test new materials inside an arc-heater based environments simulator. 3 The 
computational methodology development is achieved by implementing a solid heat transfer analysis onto an existing 
computational fluid dynamics (CFD) methodology - an Unstructured-grid Navier-Stokes Internal-external CFD 
Code (UNIC), and by simulating the heat transfer characteristics of a cylindrical pellet made of materials used in 
legacy tests that experienced difficulties, 2 heated by impinging hot hydrogen jet inside the arc-heater environments 
tester. Numerical procedures pertinent to conjugate heat transfer and physical properties germane to those legacy 
materials were implemented. The conjugate heat transfer implementation was anchored with a standard solid heat 
transfer code. The simulations of the heat transfer characteristics of the legacy materials were conducted for steady- 
state, impulsive start, and a 10 s linear ramp start conditions. 


II. Computational Methodology 

A. Computational Fluid Dynamics 

The CFD methodology is based on a multi-dimensional, finite-volume, viscous, chemically reacting, 
unstructured grid, and pressure-based formulation. Time- varying transport equations of continuity, species 
continuity, momentum, total enthalpy, turbulent kinetic energy, and turbulent kinetic energy dissipation were solved 
using a time-marching sub-iteration scheme and are written as: 
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A predictor and corrector solution algorithm was employed to provide coupling of the governing equations. A 
second-order central-difference scheme was employed to discretize the diffusion fluxes and source terms. For the 
convective terms, a second-order upwind total variation diminishing difference scheme was used. To enhance the 
temporal accuracy, a second-order backward difference scheme was employed to discretize the temporal terms. 
Details of the numerical algorithm can be found in Refs 4-8. 

An extended k-e turbulence model 9 was used to describe the turbulence. A modified wall function approach was 
employed to provide wall boundary layer solutions that are less sensitive to the near-wall grid spacing. 
Consequently, the model has combined the advantages of both the integrated-to-the-wall approach and the 
conventional law-of-the-wall approach by incorporating a complete velocity profile and a universal temperature 
profile 7 . A 2-species, 1-reaction detailed mechanism 10 was used to describe the hydrogen dissociation and 
recombination chemical kinetics. 

B. Computational Conjugate Heat Transfer in Solids * 

The solid heat conduction equation is solved with the gas-side heat flux distributions as its boundary conditions. 
The solid heat conduction equation can be written as: 


dpC p T 

dt 
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where Q v and Q s represent source terms from volumetric and boundary contributions, respectively. The present 
conjugate heat transfer model solves the heat conduction equation for the solid blocks separately from the fluid 
equations. The fluid/solid interface temperature, which is stored at interior boundary points, is calculated using heat 
flux continuity condition. For solution stability and consistency, the fluid/solid interface boundary temperature is 
updated using the transient heat conduction equation (7). The heat conduction equation on the fluid-solid interface 
can be discretized as: 
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where Q s is the wall heat flux from the solid wall to the fluid, which has an opposite sign to that from the fluid to 
the solid. The Vi factor on the left hand side of the above equation comes from the fact that only half of the solid cell 
is involved in the control volume. The following expression is obtained for updating the interface temperature. 

7o'=( 1-A)T 0 + A7j ~BQ S (9) 


where 
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Since Eq. (8) is in transient form, it is applicable to 
transient applications. For steady-state solutions, 
acceleration factor could be applied to force faster 
convergence of heat conduction in the solid. In this 
work, a multiplication factor of 10 is applied to A and 
B for steady-state solutions. 


III. Test Fixture Description 

The cylindrical specimen to be analyzed is inserted 
in a sample holder and placed in a test fixture designed 
for new materials development and testing 3 in a 
separate effort, as shown in Fig. 1. The test fixture is 
directly mated to an arc-heater (not shown) that provides 
hot-hydrogen flow. Optical ports are fitted to allow real- 
time pyrometer and laser diagnostics measurements for 
material surface temperature and centerline gas 
temperature, respectively. The copper test chamber is 
water-cooled. A tungsten/rhenium shield which serves 
as a convection shield and a radiation shield, allows the 
hot hydrogen jet to impinge on the specimen with 



Fig 1. Test apparatus. 


maximum energy and minimum heat loss to the copper chamber. UNIC code was used to design the test fixture and 
details of the design analyses can be found in Ref. 3. In this work, the computational domain includes all the flow area as 
depicted in Figs 2 and 4, except for the arc-heater and the cooling channels. The effect of the cooling channels was 
represented as wall boundary condition with a separate calculation. 3 


IV. Computational Grid Generation 

Hybrid computational grids were generated using a 
software package GRIDGEN. 11 Axisymmetric grids 
were used assuming the flow is symmetric about the 
axial axis. A series of grid verification studies using grid 
sizes ranging from 10,024, to 19,448 were performed to 
determine the current grid size (16,880). 3 It was found 
during the grid study that the computed gas-solid 
interface temperature is most favorable when structured- 
grid layers are present on both sides of the interface. 
Figure 2 shows such a grid layout where structured 
(rectangular) cells are used for the specimen and the 
region surrounding it. Structured cells are also used near 
wall boundaries, and are used for flow inlet and outlet 
boundaries, while the rest of the interior region is filled 





Fig. 2 Hybrid computational grid layout colored by 
typical computed temperature contours of hot 
hydrogen jet impinging on a rod-shaped specimen. 


with unstructured (triangular) cells (mostly not shown). It can be seen from Fig. 2 that heat is being transferred from 
the impinging hot hydrogen jet to the specimen, and the computed specimen temperature contours exhibit an 
approximately one-dimensional behavior. 
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V. Boundary Conditions, thermal conductivities and run matrix 

No-slip condition was applied to the solid walls. Fixed 
mass flow rate boundary condition was used at the inlet, 
and mass conservation boundary condition was used at the 
exit. For a conservative analysis, the hot hydrogen enters 
the chamber at 3500 K. A fixed temperature of 400 K was 
estimated for the chamber wall and the wall temperature of 
the shield facing the chamber wall was estimated to be 
2600 K, through separate one-dimensional heat transfer 
calculations. Adiabatic condition was applied to walls of 
opposing-side shield and flanges, and graphite insert (the 
convergent section). The hydrogen inlet mass flow rate 
was 10 g/s and the chamber pressure was 35 atm. A series 
of pre-calculations were performed to iterate the inlet 
temperature and species concentrations such that the inlet 
species concentrations correspond to a state of temperature 
at 3500 K for a nominal operating condition. 

Table I shows the run matrix. Benchmark cases 1, 2, and 
3 are steady-state runs with thermal conductivities of 50, 

80, and 140, respectively. The rest of the cases use variable 
thermal conductivities from measurement. Simulation 
cases 1, 2, and 3 are steady-state runs, while simulation 
cases 4, 5, and 6 are impulsive start runs, and simulation 
case 7 is a 10 s ramp start run. Composite A and B, and irradiated composite B are the legacy materials tested in the 
1970s. The irradiated Composite B is composite B material after test. The thermal conductivities of these specimens as 
function of temperatures 2 are shown in Fig. 3. Composite B represents about the highest low temperature thermal 
conductivity that can be achieved that still retains desirable mechanical properties. The thermal conductivities of the 
irradiated composite B appear quite different from those of un-irradiated composite B, and the reason of which was 
speculated as a result of annealing that occurred during the time the measurements were being taken. 2 


Table I. Run matrix. 


Case 

state 

K 

Benchmark 1 

Steady 

50 

Benchmark 2 

Steady 

80 

Benchmark 3 

Steady 

140 

Simulation 1 

Steady 

Composite A 

Simulation 2 

Steady 

Composite B 

Simulation 3 

Steady 

Irradiated composite B 

Simulation 4 

Impulsive start 

Composite A 

Simulation 5 

Impulsive start 

Composite B 

Simulation 6 

Impulsive start 

Irradiated composite B 

Simulation 7 

10 s ramp start 

Composite B 



Fig. 3 Variable thermal conductivities for three 
specimens. 


VI. Results and Discussion 
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Figure 4 shows typical computed temperature, H (atom) 
concentration contours and streamlines under nominal 
steady-state test conditions, for the entire computational 
domain. It can be seen from the temperature contours that 
the entire thrust chamber is almost uniformly heated, 
except for the region between the shield and the cooled 
chamber wall. This is because the shield serves as both 
convection and radiation shields. Nevertheless, there is still 
some heat loss to the environment and the hot hydrogen jet 
temperature drops slightly as it impinges and wraps around 
the material specimen. When temperature decreases, H 
recombines to become H 2 . This is shown in the H contours 
where its concentration decreases as the hot jet approaches 
the sample. The computed streamlines show an expanding 
hot-hydrogen jet impinging on and flowing around the 
specimen, and later exhausting into the convergent exit 
section. A large recirculation zone appears in the divergent 
section of the chamber, while a small recirculation region 
forms behind the sample holder; both of which are strongly 
affected by turbulence. Note although the plots show a 
distinctively strong hot-hydrogen jet and two recirculated 
flow regions, the pressure and Mach number contours are 
fairly uniform inside the chamber (not shown), due to the 
largely low subsonic flow field and the protection of the 

shield. The flow becomes supersonic in the exit pipe. Fi S- 4 From l »P botl « m ’ ^ computed steady-state 

temperature and H Mass Fraction Contours and Streamlines. 

Steady-State Benchmark Analysis with Constant Thermal Properties 

Figure 5 shows a comparison of the computed steady-state centerline coupon temperatures, using conjugate heat 
transfer methodology developed in this work, 
with those obtained from a standard solid heat 
transfer code SINDA . 12 Solutions for three 
constant thermal conductivities, 50, 80, and 140, 
measured at room temperature for flow element 
materials heat treated at 2800, 3050, and 3 120 K, 
respectively , 2 were performed. The first iteration 
of SINDA solution was based on three arbitrarily 
chosen wall temperatures of 300, 2500, and 3400 * 

K, by assuming the final solid temperatures are 
between 300 and 3400 K. It should be noted that 
these fixed wall temperatures imply three wall 
heat flux profiles. It can be seen that the UNIC 
computed centerline temperatures compared 
reasonably well with those obtained from the first 
SINDA iteration, and the maximum deviation is 
less than 10 K. Based on the temperatures 
calculated in the first iteration, the second SINDA 
iteration was performed and the three input wall 
temperatures were narrowed to 3405, 3410, and 
3413 K. It can be seen that solutions of UNIC and pjg 5 a comparison of computed steady-state centerline 

those of second SINDA iteration are now very spec imen temperatures, 

well compared, and the deviation is less than 2 K. 

In summary. Fig. 5 demonstrates that the UNIC computed solid temperatures compare very well with those of a 
standard solid heat transfer code SINDA. Both UNIC and SINDA solutions exhibit similar behavior. For example, 
high thermal gradients occur near the stagnation point (x = 0). In addition, the lowest thermal conductivity (k = 50) 
gives the highest thermal gradient, while the highest thermal conductivity (K = 140) results in the lowest thermal 
gradient. 
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Analysis with Variable Thermal Properties 


Steady-State Analysis 

Three constant thermal conductivities were used to 
perform benchmark analysis as mentioned above. For 
actual flow element materials such as Composite A, 

Composite B and irradiated Composite B, their 
conductivities varies with temperature, as shown in Fig. 

3. Figure 6 shows the computed steady-state centerline 
temperatures for these three specimens, simulated as 
cylindrical pellets heated by hot hydrogen jet 
impingement as described in Fig. 1. It can be seen that 
the centerline temperature profiles for Composite A and 
Composite B overlap, while that for the irradiated 
Composite B deviate only slightly near the stagnation 
point (x = 0) and near x = 0.04. This is because the 
cylindrical pellet is surrounded by hot hydrogen flow 
and the whole pellet is heated up to temperatures from 
about 3400 to 3425 K, as shown in Figs. 2 and 4. At 
these temperature levels, the thermal conductivities are 
nearly constant (30 ~ 35 W/mK), as shown in Fig. 3. 

The temperature scale in Fig. 6 is zoomed in, that is, Fig. 6 Computed steady-state centerline temperatures 
the steady-state thermal gradients for all three specimens f° r three hypothetical specimens, 
appear to be benign, implying that there was no high 

thermal gradient or high thermal stress for these pellets, during the steady-state heating by hot hydrogen. Although the 
set up of impingement heating of a short pellet (fig. 1) is not intended to explain the mid-section thermal anomaly of a 
internally heated, full-length flow element, rather a simpler way to isolate the physics. With the range of thermal 
conductivities shown in Fig. 3, it is speculated however, that there should be stronger thermal gradients during the 
transient impingement heating process. That possibility is investigated by performing three impulsive and one linear 10 s 
ramp startup computations, as discussed in the following. 

Transient Analysis - Impulsive Start 

Three impulsively started transient conjugate heat 
transfer computations were performed for the three 
specimens with the variable thermal conductivities as 
described in Fig. 3. Hydrogen gas was assumed to be in 
the tester initially at initial temperature and pressure of 
300 K and 1 atm, respectively. The initial pellet * 
temperature was assumed to be at 300 K. For the most ^ 
conservative case, hot hydrogen was assumed to enter 
the tester impulsively at 3500 K and 35 atm, 
respectively, meaning the inlet hydrogen gas was heated 
to 3500 K and pressurized to 35 atm instantly by the arc 
heater at time zero. 

As mentioned in the last section, it was anticipated 
that impulsive would create high thermal gradient, or 
thermal shock in the solids. Figure 7 shows the 
computed transient centerline temperature profiles for 
composite A. It can be seen that unlike the steady-state FJg ? Computed cy i indrical specimen centerline 

profile, strong thermal gradients occur at both ends of temperature history for Composite A during an 

the pellet, where the impingement end or the stagnation i mpu lsive start 
point end (x = 0) encounters the strongest thermal 
gradient throughout most of the transient heating 

process, especially during the 3 to 7 s period. During that time period, it can be seen that the center core had a difficult 
time to get heated up from 1600 to 1800 K, resulting in strong thermal gradients at the two ends. This phenomenon is 

7 

American Institute of Aeronautics and Astronautics 




x, m 





strongly correlated to the low thermal conductivity at high temperatures for composite A, as shown in Fig. 3. 

Figures 8 and 9 show the computed transient temperature profiles for composite B and irradiated composite B, 
respectively. It can be seen that except in the early heating period (0.05 ~ 1.5s) when the irradiated composite B shows 
stronger thermal gradient than that of composite B, both show similarly strong thermal gradients at the two ends when 



x, m 


x, m 


Fig. 8 Computed cylindrical specimen centerline Fig. 9 Computed cylindrical specimen centerline 

temperature history for Composite B during an temperature history for irradiated Composite B during 

impulsive start. an impulsive start. 

the central core reaches 1600 K. In fact, during that time period, the transient temperature profiles of composite B and 
irradiated composite B look similar to those of composite A (Fig. 7). k 

The similarity in the occurrence of high thermal gradient during the 3 ~ 7 s time period fori composite A and B, and 
irradiated composite B correlates to the fact that the thermal conductivity of all three specimen converge to a similarly 
low value as temperature exceeds 1,600 K, even though their individual thermal conductivity profiles look quite different 
in lower temperatures. Figures 7, 8, and 9 demonstrate that the thermal wave propagation inside the cylindrical pellet is 
a strong function of the characteristics of its thermal conductivity. In addition, it could be inferred that for other 


specimens with completely different thermal conductivity 
characteristics, different thermal transient histories would 
have been computed. Figures 7, 8, and 9 also show that, 
for the three cylindrical flow element specimens studied, it 
takes more than 12.5 s to heat the pellets up, using a 
impulsive start. 

Transient Analysis - 10 s Ramp Start 
The 10 s ramp start assumes the hydrogen gas inlet 
temperature and pressure are brought linearly to 3,500 K 
and 35 atm, from 300 K and 1 atm, respectively, during a 
10 s period. This means that the leading edge of the 
cylindrical pellet does not see the final temperature until 
the 10 th s, comparing to sensing the final temperature at the 
start command of the impulsive start. Only the effect of 
thermal conductivity characteristics of Composite B was 
investigated, since the impulsive start analyses show that 
the transient behaviors of all three specimens are similar at 
high temperatures. Figure 10 shows the computed transient 
centerline temperature history for Composite B during the 
10 s ramp start. It can be seen that from 0.05 to 8.5 s, the 
thermal gradients are fairly benign at both ends of the 



0 0.01 0.02 0.03 0.04 0.05 0.06 

x, m 


Fig. 10 Computed cylindrical specimen centerline 
temperature history for composite B during a 10 s 
ramp start. 
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cylindrical pellet. At 9.5 s, as the center core temperature reaches 1,600 K, higher thermal gradients again start to show. 
Hence, the thermal conductivity effect from Fig. 3 still dominates in the 10 s ramp start, even thought the high thermal 
gradient occurs about 6.5 s later than those of the impulsive starts. The total heat up time is more than 15.5 s. 


VII. Conclusion 

A computational conjugate heat transfer methodology was developed, to study the heat transfer effect of a hot 
hydrogen jet impinging on a cylindrical pellet made of hypothetical thermal nuclear rocket flow element materials, 
as a first step towards developing a multiple physics, thermo-fluid computational methodology to predict 
environments for hypothetical solid-core, nuclear thermal engine thrust chamber and components. This is 
accomplished by implementing a solid heat transfer procedure onto an existing computational fluid dynamics 
formulation. The solid heat transfer module was anchored and the computed pellet centerline temperatures compare 
very well with those of a standard solid heat transfer code. Analyses were then performed for hypothetical 
cylindrical pellet using variable thermal conductivities of three legacy flow element materials. Under the 
assumptions made in this study, it was found that high thermal gradients occur during the transient heating processes 
and not the steady-state heating processes. High thermal gradients occur at both ends of the pellet, especially near 
the impingement point. The transient heating behavior inside these rod-shaped specimens depends solely on their 
thermal conductivity characteristics. For the three specimens studied, since the thermal conductivities of all three 
specimen converge to a constant low value as temperatures exceed 1,600 K, high thermal gradient starts to form at 
both ends of the pellet when the center core temperatures reaches 1,600 K, during either the impulsive or 10 s ramp 
start process. 
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The objective of the Multi-physics Thrust Chamber Modeling task is 
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Introduction - continued 
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Multi-physics invoked include hydrogen dissociation kinetics and 
thermodynamics, turbulent flow, and conjugate heat transfer. 


UNIC CFD code 







UNIC Computational Methodology 
Conjugate Heat Transfer: Gas-side Governing Equations 



Gas-side Governing Equations - Continued 





Conjugate Heat Transfer: Solid-side Governing Equations 
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Hot Hydrogen Materials and 
Component Tester 



Optical Port 



NASA MSFC Multi-Gas Arc-Heater Facility 
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nuclear thermal environment simulator 



Hypothetical 


*— Tl 

g C 
§ ° 

JL O) 

3 -j= 

LL C3 


£ ^ CO 

^§o 

LO LO ° 
CO CO ^ 


cl h E 


O O rf r R -Q H 
cri OO i-H U U 


<J pq O <tj PQ 

CD CD <D CD 

mm 

# |“bi 4 *i> jssfsl 

co co CD 
O O tg 
Cl pL «2 

■ S3 . S HD 

8 ' M- g 

O O ?H 




T 3 

>> 

kQ 

kQ 

K*> 


t 3 

03 

*“d 

<5 

CD 

CO 

CD 

,j ja ^ 

§ 

cO 

CD 

cc 3 

CD 

So 

m 

So 

GO 

DO 



***4 (N CO | 
^j j 

fel U S o 


CD O Q 


cN co 


9 cd c3 c<d 

“5 r ~j *~~i r ~~i 

O 3 3 r- ; 


u pci pci n « 


GO GO 100100 














Computational Grid Layout and Computed Temperature 
Contours around the Cylindrical Specimen 


mm 





m 


I 




mm 


S J Basis BH Silfil III Mil I9i{ till llli til! Mil Hill 1 1 3 1 1 1 S 1 1 S 8 S B S 8 s 



ii 

ii 

m 




■”d u a e a a 

RBHB3 

sssss 



i a ^>r a ass si ismnni mi mi mi mi mi mi mis iisiiihi iiHHBB^r^S-^abl 

L-j; iBIIIilil 3121 ilti (HI HU IHI HH Hill IISHIRl 1 8 0 1 1 ■ few ■ - 

if r** t<- ±- IBMailll till lilt llll llll EIII lilt IIEIt 11311111 1 1 9 ■ B ■ 

lanimi iiii mi mi mi mi mi him iisiiiii iisiii^st^&iy; 
w feffilSBiBiBi iBtiiiin iiii mi iiii ill! iiii im iiiii 11111111 iiiiiigBrfj^S™ 
J8I3IBI HUSIUIIIIIIIIUIIIWHIIIIIIIIIIIIIIIIIIII IIIIB1 gf 
iianesi ibiiiiiii 1111 1111 mi ill! 1111 iiii mu iianiBi iiimbi> 
ieillil 18IIBIII1 llll ill! til! IIII !IEI (ID IIIII IIIIIIII IIIIBI 
J 1 SIBI 8 I IBIIIIIIS fill HI! nil IIIIIIII llll IIIII 1IBIIIII IIIIIlL 
J&8M1I IBIIIIIII llll llll llli llll EIII llll iilllllBIIIII 118118*1 
JpIllBI llllllllllillliil llll llli IIIIIIII llllllllillll 111111%; 
llllllllillll mi llll llll llll llll Hill IIIIIIII IIHSlf 
JV-*S»££i 8 fiBJSii iiitimi nil mi mi nil mi mi iiiii iibiiiii usisia 
4Bisiii iiiiiiiiuiiiimimiMiiiuiiiuimuiiiin iibibi* 

.. • ijipijiiiiiiiijiiiiiii mini 

iMiiBiilifimiJiCiiaifial ns 1 a i ZS£ 

piiisir^iiiSimiiiaiiiil wiibiST 



§ * p- ^tl a b sisti Jiaiif rnmm m 


pamam 

i*0ipasi nmmm 
kibb sa^i mumsm 

anstaxt irastxjiyt,' 

mmnatm 

amaxiwi 






. 





?#38BD*SJ: 
!’*en k b •• ss ' 
v*V 



111 


3 &L 

m 

fli 

ii 

SK» Hi 

fc I 

sg 

B 
ii 

SSb 




mHnwnMnMfnnmi 
if 1 ■ 1 iBsisiaillifilii 

•‘tin ihj-» 1 1 ■ 1 lifl’-ii in mill 1 bib 1 gswgs 

w«a:- .*^r w^ii^piffiniiiiiiisSiK 

1 . 1 .li'iivjii nniiiiusriiiipMMismjiiiiiiii^^aS 

# M lUmummmXmmiliimffsm 


I m 

j&tF 


MlllHII^ 

lifipiiiiiiigiSiK. 
Bliliili^pl! 


Silfiiilsff 

$ 8(11911 
*ili! 81 l 
IIIIIIII 
*11313111 
IIIIIIII 
HilllSIE 

MHIllliSfllfPf 

isiHiiiiiIiSiiiii; 




jhbb bih-hh 
iififiniTi^iiHiiiiiiQiL, 
liiBli iiii 1111 1111 ill 
liiffipliiiiiiii iiii iii 
MH tp nil iiiuiii iiiii 

rffiiiii! 

iiiWiiifi 1111 1111 mi ini iinpiiuiits 1 1 a i a i 
Jim mil I!!! III! ill! llliiMiJIlillllll 
iiEliiiisiii! 

jflljlllllli 

IIIIIUHUIIlBllWillllllli 



till llll llll llli llll ill) l!tl|BMlllll? 3 tf^ 




nl 



Typical Steady-State Computational Results 




Hot Hydrogen Materials and Component Tester 
Conjugate Heat Transfer Analysis - UNIC & SINDA Comparison 
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SINDA solutions 








Computed cylindrical specimen centerline temperature history 
for composite A during an impulsive start 
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Animation of transient temperature contours of a cylindrical composite A 
heated by hot hydrogen using an impulsive start 






Computed cylindrical specimen -centerline temperature history 
for composite B during ah impulsive start 
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Computed cylindrical specimen centerline temperature history 
for irradiated composite B during an impulsive start 
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Computed cylindrical specimen centerline temperature history 
for composite B during a 10 s ramp start 
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The transient heating behavior Inside these rod-shaped specimens 
depends solely on their thermal conductivity characteristics. 


